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Abstract 

We discuss the dynamics of a Bose-Einstein condensate in a double- well trap subject to phase 
noise and particle loss. The phase coherence of a weakly-interacting condensate as well as the 
response to an external driving show a pronounced stochastic resonance effect: Both quantities 
become maximal for a finite value of the dissipation rate matching the intrinsic time scales of the 
system. Even stronger effects are observed when dissipation acts in concurrence with strong inter- 
particle interactions, restoring the purity of the condensate almost completely and increasing the 
phase coherence significantly. 
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In our naive understanding thermal noise is generally deconstructive, hindering measure- 
ments and degrading coherences in quantum mechanics. A paradigmatic counterexample to 
this assertion is the effect of stochastic resonance (SR), where the response of a system to 
an external driving assumes its maximum in the presence of a finite amount of thermal noise 
This maximum occurs when the time scales of the noise and the driving match. By 
now, SR has been shown in a variety of systems, an overview is given in the review articles 

In addition to numerous examples in classical dynamics, SR has also been found in a 
variety of quantum systems (see 5] and references therein). Recently, there has been an 
increased interest in controlling and even exploiting dissipation in interacting many-body 
quantum systems. For instance, the entanglement in a spin chain assumes an SR-like max- 
imum for a finite amount of thermal noise 6j. Furthermore, it has been shown that dissi- 
pative processes can be tailored to prepare arbitrary pure states for quantum computation 
and strong! 



turn gates 



correlated states of ultracold atoms [7] or to implement a universal set of quan- 
. Actually, a recent experiment has even proven that strong inelastic collisions 
may inhibit particle losses and induce strong correlations in a quasi one-dimensional gas of 
ultracold atoms j^]. 

In this letter we demonstrate the constructive effects of dissipation such as stochastic 
resonance for an interacting many-particle quantum system realized by ultracold atoms 
in a double-well trap with biased particle dissipation. It is shown that a proper amount of 
dissipation increases the coherence of the atomic cloud, especially in concurrence with strong 
inter-particle interactions. These effects are of considerable strength for realistic parameters 
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and thus should be observable in ongoing experiments [10, llll, LUJ H3[ • 

The unitary dynamics of ultracold atoms in a double- well trap is described by the two- 
mode Bose-Hubbard Hamiltonian 
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where Oj and aj are the bosonic annihilation and creation operators in the jth well and 
hj = a,jdj are the number operators. In general we consider a biased double-well trap, where 
the ground state energies of the two wells differ by 2e. We set h = 1, thus measuring all 
energies in frequency units. 



The main source of decoherence is phase noise due to elastic collisions with atoms in 



16] which effectively heats the system. The heating rate is fixed 



the thermal cloud [15|, 

as 7 g = 5s _1 in the following, which is a realistic value for the experiments in Heidelberg 



10, Methods to attenuate this source of decoherence were discussed only recently 

17| . Amplitude noise, i.e. the exchange of particles with the thermal cloud due to inelastic 
scattering, drives the system to thermal equilibrium. However, this effect is usually much 
;oo weak to produce the effects discussed below in present experiments (cf. the discussion in 
Jil). In contrast, a strong and tunable source of dissipation can be implemented artificially 
by shining a resonant laser beam onto the trap, that removes atoms with the site-dependent 
rates 7 a j from the two wells j = 1,2. Non-trivial effects of dissipation such as the stochastic 
resonance discussed below require strongly biased loss rates, i.e. j a i ^ 7 a2 . For a laser beam 
focused on one of the wells an asymmetry of f a = (j a 2 - 7oi)/(To2 + 7oi) = 0.5 should be 
feasible. Thus we consider the dynamics generated by the master equation 

i=i,2 

~2 Yl 7a J + P h % ~ 2a iP a l) • ( 2 ) 

The macroscopic dynamics of the atomic cloud is well described by a mean-field approxi- 
mation, considering only the expectation values Sj = 2tr(Ljp) of the angular momentum 
operators 
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and the particle number n = tr((rii + n2)p)- The time evolution of the Bloch vector s and 
the particle number is then given by 

S X = IdSy USyS Z ?2 s X) 

Sy = 2,.J S z + 1tS X ~\- U S X S Z Sy, 

S z = 2 J Sy S z fa n i 

n = ~T{ l n-T{ l f a s z . (4) 

As usual expectation values of products have been factorized in the [/-dependent interaction 
terms to obtain a closed set of evolution equations whereas the dissipation terms 



FIG. 1: (Color online) Contrast a in the quasi-steady state in dependence on the tunneling rate J 
and the dissipation rate 1/T\ for U = and e = 0. 



are exact. Furthermore we have defined the transversal and longitudinal damping times by 

Tf 1 = ( 7 ai+7a 2 )/2 and T^ = lp + T{\ (5) 

These equations of motion resemble the celebrated Bloch equations in nuclear magnetic 
resonance with some subtle but nevertheless important differences. The longitu- 

dinal relaxation is now associated with particle loss and, more important, the dynamics is 
substantially altered by the interaction term [lo|, HQ. 

In the following we will show that a finite amount of dissipation induces a maximum 
of the coherence which can be understood as an SR effect. In this discussion we have to 
distinguish between two different kinds of coherence, which will both be considered in the 
following. First of all we consider the phase coherence between the two wells, which is 
measured by the average contrast in interference experiments as described in flpl . \\\ and 
given by 

® = 77 t~t = • (o) 

(ni + n 2 ) n 

Secondly, we will analyze how close the many-particle quantum state is to a pure Bose- 



Einstein condensate (BEC), which is a coherent state for the SU(2) operator algebra 22]. 
This property is quantified by the purity p = 2tr(p 2 cd ) — 1 = |s| 2 /n 2 of the reduced single- 
particle density matrix p re d, cf. [20 ] . 

Let us first discuss the weakly-interacting case, where the mean-field equations of motion 
(j4|) provide an excellent description of the dynamics, which is exact for U — 0. Obviously, 
only the trivial solution s = and n = is a steady state in the strict sense. However, 
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FIG. 2: (a) Average contrast a after 1.5 s of propagation starting from a pure BEC (i.e. a product 
state) with s z = re/2 and re(0) = 100 particles in dependence of the tunneling rate J for T\ = 0.5 s, 
e = 10 s -1 , U = 0.1 s -1 . (b) Histogram of the probabilities to measure the relative phase (j) and 
the relative population imbalance s z /n in a single experimental run after t = 1.5 s obtained from a 
MCWF simulation of the many-body dynamics. 

the system rapidly relaxes to a quasi-steady state where the internal dynamics is completely 
frozen out and all components of the Bloch vector and the particle number decay at the 
same rate. Fig. [1] shows the contrast a for this quasi-steady state as a function of the 
tunneling rate J and the dissipation rate 1/Xi for U = 0. For a fixed value of one of the 
parameters, say J, one observes a typical SR-like maximum of the contrast for a finite value 
of the dissipation rate 1/Ti. In particular, the contrast is maximal if the time scales of the 
tunneling and the dissipation are matched according to 24] 



f a T{ 1 « 2 J. 



(7) 



This scenario is robust and not altered by weak inter-particle interactions. Changes in the 
system parameters such as e preserve the general shape of a (1/Ti, J) and the existence of 
a pronounced SR-like maximum. At the most, the function a (1/Ti, J) is stretched, shifting 
the position of the SR-like maximum. 

The occurrence of a maximum of the contrast is explained by Fig. [2] (b), where the results 
of a Monte Carlo wave function (MCWF) simulation 23( of the many-body dynamics are 
shown for three different values of J and U = 0.1s -1 . We have plotted a histogram of 
the probabilities to observe the relative population imbalance s z and the relative phase </> 




FIG. 3: (a) Oscillation of the relative population imbalance s z /n of a weakly driven two-mode 
BEC for Jo = 1.5 s , T\ = 0.5 s and e = 0. (b) Amplitude of the oscillations in dependence on the 
tunneling rate Jo- 

in a single experimental run for three different values of the tunneling rate J after the 
system has relaxed to the quasi-steady state. With increasing J, the atoms are distributed 
more equally between the two wells so that the single shot contrast increases. Within the 
mean-field description this is reflected by an increase of v/s^~+s^/|s| at the expense of \s z \. 
However, this effect also makes the system more vulnerable to phase noise so that the relative 
phase of the two modes becomes more and more random and |s|/n decreases. The average 
contrast ([6]) then assumes a maximum for intermediate values of J as shown in Fig. [2] (a). 
Note that the trap is assumed to be weakly biased in this example, shifting the position 
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So far we have demonstrated a stochastic resonance of the contrast for a BEC in a 
static double-well trap with biased particle losses. In the following we will show that the 
system's response to a weak external driving also assumes a maximum for a finite dissipation 
rate - an effect which is conceptually closer to the common interpretation of stochastic 
resonance. We consider a weak driving of the tunneling rate J{t) = J + J\ cos(ut) at the 
resonance frequency u> = \/ Jq + e 2 , where the amplitude is not more than Ji/Jo = 10%. 
This can be readily implemented in optical setups by varying the intensity of the counter- 
propagating lasers forming the optical lattice. Fig. [3] (a) shows the resulting dynamics 
for Ti = 0.5 s and Jq = 1.5 s" 1 . After a short transient period, the relative population 
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FIG. 4: Time evolution of the purity p and the contrast a for J = U = 10 s -1 , e = 0, T\ = 0.5 s. 
The initial state is a pure BEC with s z = n/2 and n(0) = 100 particles. The results of a MCWF 
simulation averaged over 100 runs are plotted as a thin solid line while the mean-field results are 
plotted as a thick line. The dashed line shows the steady state values for 1/T\ = I/T2 = 0, i.e. 
without coupling to the environment. 

imbalance s z (t)/n(t) oscillates approximately sinusoidally. The system response measured by 
the amplitude of these forced oscillations shows the familiar SR-like maximum as illustrated 
in Fig. [3] (b). It should be detectable without major problems in ongoing experiments, in 



which the population imbalance s z can be measured with a resolution of a few atoms 



3, 



A more detailed study of the parameter ranges for such a driven case of SR will be discussed 



in a forthcoming article 



24]. 



Even more remarkable values of the coherences are observed in the case of strong interac- 
tions, which is experimentally most relevant and theoretically most profound. The interplay 
between interactions and dissipation significantly increases the coherences in comparison 
to situations where one of the two is weak or missing. An example for the dynamics of a 
strongly-interacting BEC is shown in Fig. H]for an initially pure BEC with s z = n/2, calcu- 
lated both with the MCWF method and within the mean-field approximation At first, 
the purity p and the contrast a drop rapidly due to the phase noise and the interactions, 



cf. 



20| | . For intermediate times, however, the system relaxes to a nonlinear quasi-steady 



state, which is a nearly pure BEC mostly localized in the well with the smaller decay rate. 
Consequently, the purity p is restored almost completely and the contrast a is relatively 
large. In close analogy to the celebrated self-trapping effect 0, Q, 19], this quasi-steady 
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FIG. 5: Purity p and contrast a after t = 2 s in dependence on the dissipation rate \jT\ calculated 
within the mean-field approximation. The remaining parameters are chosen as in Fig. [5J 



state exists onl y a s long as the effective interaction strength Unit) is larger than a critical 
value given by [25j 

U 2 n 2 > A J 2 - f a T{ 2 . (8) 

As the particle number n decays, this state ceases to exist so that the system relaxes to a 
linear quasi-steady state with considerably smaller values of p and a as discussed above for 
the weakly-interacting case. 

Moreover, the coherences at intermediate times are also larger than in an interacting, 
but non-dissipative system. The dashed lines in Fig. [4] show the steady state values of the 
purity p and the contrast a for 1/Ti = 1/T 2 = 0, apart from occasional revivals due to the 
finite particle number. It is observed that the coherences are considerably smaller compared 
to the strongly-interacting open system. This loss of coherence can be understood by the 
fact that the interactions lead to an effective decoherence on the single-particle level 20] . 
degrading a and p. This effect is mostly cured by the dissipation. 

The behaviour illustrated in Fig. H] and discussed above is universal in the sense that 
the maxima of the purity and the contrast are present for all values of U and 1/Ti if only 
Un(t = 0) is well above the critical value ([Sj) for the existence of the nonlinear quasi-steady 
state. However, the maxima occur later if Ti or U increase. The purity p and the contrast 
a after a fixed time t = 2 s are plotted in Fig. [5] in dependence on the dissipation rate 1/Ti, 
showing pronounced maxima for finite values of 1/Ti. For smaller dissipation rates, the 
maximum of the contrast has not been assumed yet while the system has already relaxed to 
the linear quasi-steady state for larger values of 1/Ti. 

To summarize, we have shown that the coherence properties of a weakly and in partic- 
ular also of a strongly interacting Bose-Einstein condensate in a double-well trap can be 
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controlled by engineering the system's parameters and dissipation simultaneously. An im- 
portant conclusion is that the interplay of interactions and dissipation can drive the system 
to a state of maximum coherence, while both processes alone usually lead to a loss of coher- 



ence. Since the double-well BEC is nowadays routinely rea 



ized with nearly perfect control 



on atom-atom interactions and external potentials [10|, llll], we hope for an experimental 
verification and future extensions of the predicted stochastic resonance scheme. 
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